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We study models of a society composed of a mixture of conformist and reasonable contrarian 
agents that at any instant hold one of two opinions. Conformists tend to agree with the average 
opinion of their neighbors and reasonable contrarians to disagree, but revert to a conformist behavior 
in the presence of an overwhelming majority, in line with psychological experiments. The model 
is studied in the mean field approximation and on small-world and scale-free networks. In the 
mean held approximation, a large fraction of conformists triggers a polarization of the opinions, 
a pitchfork bifurcation, while a majority of reasonable contrarians leads to coherent oscillations, 
with an alternation of period-doubling and pitchfork bifurcations up to chaos. Similar scenarios are 
obtained by changing the fraction of long-range rewiring and the parameter of scale-free networks 
related to the average connectivity. 
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I. INTRODUCTION 

In a previous work [I], we studied the collective be¬ 
havior of a society of reasonable contrarian agents. The 
rationale was that in some cases, and in particular in 
the presence of frustrated situations like in minority 
games [2], it is not convenient to always follow the major¬ 
ity, since in this case one is always on the “loosing side” of 
the market. This is one of the main reasons for the emer¬ 
gence of a contrarian attitude. On the other hand, if all 
or almost all agents in a market take the same decision, 
it is often wise to follow such a trend. We can denote 
such a situation with the word “social norm”. Following 
an overwhelming majority is an ecological strategy since 
it is probable that this coherent behavior is due to some 
unknown piece of information, and in any case the com¬ 
petitive loss is minimal since it equally affects the other 
agents. Indeed, it is well known that in the presence of 
an overwhelming majority, individuals tend to align with 
it, even if it is in contrast to evidence [3j. 

The agents in our model can hold opinion 0 or opinion 
1, representing one of two parties, or against/in favor of 
an option. Each agent gathers the average opinion of 
his neighbors and changes his opinion according to this 
average. The time evolution is synchronous, so the model 
is essentially a cellular automaton. 

Contrarians were introduced in a different socio¬ 
physical model by S. Galam HE]. In this case they 
had the effect of destroying consensus in a society mainly 
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formed by conformists. In this case, when the fraction of 
contrarians becomes opinion-dependent, chaotic dynam¬ 
ics appear [6]. 

A society fully composed by reasonable contrarians ex¬ 
hibits interesting behaviors when changing the topology 
of the connections. On a one-dinrensional regular lattice, 
there is no long-range order, the evolution is disordered 
and the average opinion is always halfway between the 
extreme values 0 and 1. However, adding long-range con¬ 
nections or rewiring existing ones, we observe the Watts- 
Strogatz “small-world” effect, with a transition towards a 
mean field behavior. But since in this case the mean field 
equation is, for a suitable choice of parameters, chaotic, 
we observe the emergence of coherent oscillations, with 
a bifurcation cascade eventually leading to a chaotic-like 
behavior of the average opinion. The small-world transi¬ 
tion is essentially a synchronization effect. Similar effects 
with a bifurcation diagram resembling that of the logistic 
map have been observed in a different model of “adapt if 
novel - drop if ubiquitous” behavior, upon changing the 
connectivity pm 

Since a homogeneous society of unreasonable contrar¬ 
ians is not so reasonable, we study here the problem of 
collective behaviors in the presence of a mixture of con¬ 
formists and contrarians. To keep things simple, rea¬ 
sonable contrarian and conformist agents have the same 
behavior in the presence of an overwhelming majority 
of their neighbors. Conformists become less conformists 
in the presence of a large majority. We can call them 
“slightly unreasonable”. The presence of reasonable con¬ 
trarians and slightly unreasonable contrarians avoids ab¬ 
sorbing states which are rather unusual in real societies. 

In the presence of a strong fraction of conformists, we 
have the classical ferromagnetic Ising scenario, with the 



appearance of a stationary average opinion different from 
one half. As this fraction becomes smaller, this polarized 
opinion vanishes, as expected. What is unexpected is 
another bifurcation, with the appearance of oscillations 
and chaos as the fraction of conformists becomes smaller. 

The outline of the present paper is as follows: in Sec.[TT| 
we present the model in detail. Its mean field approxima¬ 
tion is discussed in Sec. EH Then the model is studied 
on small-world networks, Sec. IV and on on scale-free 
networks, Sec. |V|. W e end the presentation with some 
conclusions, Sec. I VII 



h 


II. THE MODEL 


Our model society is formed by N agents with a 
fraction £ of conformists and a fraction 1 — £ of con¬ 
trarians. Agent i, i = 0,..., N — 1 holds an opinion 
s(i,t) £ {0,1} at time t. The opinions of all agents 
change synchronously in time. Agent i gathers the aver¬ 
age opinion of his neighbors and changes opinion tending 
to agree with his neighbors if he is a conformist, or to 
disagree if he is a contrarian. 

The neighborhoods of all agents are defined by the ad¬ 
jacency matrix A with components £ {0,1} in a way 
that dij = 1 if agent j belongs to V s neighborhood and 
dij = 0 if he doesn’t. The connectivity ki of agent i is 
the number of agents in i’s neighborhood, 

h = (!) 

0 

and the average opinion hi of his neighbors is 


hi = 7 -^2 a ijS{j)- 


The average opinion c of the society is 

1 




( 2 ) 


(3) 


Given the average opinion h of the neighbors of agent 
i at time t , s(i,t + 1 ) = 1 according to the transition 
probability r(h; J) defined by [T] 


! e if h < q, 

- -t—- —7 if q < h < 1 — q, 

1 + exp(- 2 J( 2 /i - 1 )) 

1 — £ \{ h > l — q. 

_( 4 ) 

In this last expression J represents the agent’s conviction, 
conformists have J > 0, reasonable contrarians, J < 0. 
The graphs of r are shown in Fig. [I] A reasonable con¬ 
formist or a contrarian can assume opinion 1 with prob¬ 
ability e (1 — e) when 0 < h < q (1 — q < h < 1). When 
q<h<l~q & reasonable conformist will probably 
agree with his neighbors, and a contrarian will proba¬ 
bly disagree. Both conformists and contrarians share the 


FIG. 1 : (Color online) The transition probabilities r given 
by Eq. © with | J| = 3 . For q < h < 1 — q, r is an increasing 
function of h (in red) for conformists, J = 3 , and a decreasing 
one (in blue), J = — 3 , for reasonable contrarians. In this and 
the folowing Figs, q = 0.1 and e = 0.2 unless stated otherwise. 


same values of | Jj, q , and e. A nonzero value of e avoids 
the presence of absorbing states but causes the “slightly 
unreasonable” behavior of conformists. Unless otherwise 
stated, we always use the value q = 0.1 and e = 0.2, as 
in Ref. P. 

In the mean field approximation c changes determinis¬ 
tically and we can characterize the properties of its tra¬ 
jectory by means of the Lyapunov exponent A. On the 
other hand, c changes probabilistically on networks, so we 
use Boltzmann’s entropy 77 piano], which is applicable 
in both cases. To define 77 we partition the unit interval 
in L disjoint equal sized subintervals A, i = 1,..., L and 
find the probability < 7 ,; that c falls in subinterval A ■ Then 

-1 L 

71 = i^z5> log(? i- ( 5 ) 

i=l 

The probabilities qi are found numerically by finding the 
frequency with which an orbit visits each subinterval A 
after a transient. A fixed point of the trajectory corre¬ 
sponds to 77 = 0 , and when the orbit visits every subin¬ 
terval A with the same frequency 77 = 1. If L = 2 b and 
the orbit is periodic with period 2 a , 77 = a/6. In the limit 
L —> 00 , 77 —>• 0 for periodic orbits. 

For deterministic maps we can compare the behavior 
of the Lyapunov exponent A and that of Boltzmann’s en¬ 
tropy 77. Values of A > 0 are equivalent to 77 > 1/2 when L 
is sufficiently large. In other words, deterministic chaos 
corresponds to 77 > 1/2 and order to 77 < 1 / 2 . By ex¬ 
tending this correspondence to the probabilistic network 
dynamics, we define disorder whenever 77 > 1/2 and order 
when 77 < 1/2. 


III. MEAN FIELD APPROXIMATION 

We start with a model of a society where the neigh¬ 
borhood of each agent i, either conformist or contrarian, 
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FIG. 2: (Color online.) Return maps of the mean field ap¬ 
proximation, Eq. §, smooth curves (in red), and 40 iterates 
of the map starting from a random initial c, stepping lines (in 
blue), (a) £ = 0.06. (b) £ = 0.2. (c) £ = 0.4. (d) £ = 0.9. 
In this and the following Figs., | J\ = 5, k = 20, e = 0.2, and 
q = 0.1 unless stated otherwise. 


is formed by k random neighbors, i.e., the mean field 
approximation for a fixed connectivity k. With a frac¬ 
tion £ of conformists and a fraction 1 — £ of reasonable 
contrarians, the time evolution of the average opinion c 
is 




FIG. 3: (Color online.) The mean field average opinion c 
as a function of the fraction of conformists £ with |J| =5, 
k = 20, and q = 0.1. (a) e = 0.2, £ c = 0.0665. (b) e = 0.0, 
£ c = 0.0795. In (a) and (b), £ a = 0.3245 and & = 0.6755 
so that 1 — £(, = £ a . For £ slightly larger than £ c the first 
and third branches (starting from the bottom) correspond to 
Co = 0.9 (in red) and the other two to co = 0.1 (in blue). For 
< £ < 1, the upper branch (in red) corresponds to Co = 0.9, 
the lower one (in blue) to Co = 0.1, For all values of £, two 
consecutive iterations are plotted after a transient of 1,000 
time steps for both Co = 0.9 (in red) and Co = 0.1 (in blue). 


Co > 1/2 and to c < 1/2 when co < 1/2 with co the 
average opinion at t = 0, (see also Fig. [2]). 

We prove in the Appendix that £ Q = 1 — £b and find 
that the approximate behavior of £ a as a function of J 
and k is given by 


£a{J,k) 



(7) 


c' = £ 


ii;=0 


’(l-c) 


k—v 


[£,r{w/k\ J) + (l-£)T(w/k;-J)], (6) 


with c and c' the average opinions at times t and t + 1 
respectively. In the right-hand side term, the first paren¬ 
thesis is the w combinations from a set of k elements. 

In Fig. [2] we show return maps of Eq. ([6]) for different 
values of £. For small £, Fig. [ 2 ] (a), the map is chaotic, 
and for larger values of £, we find periodic orbits or fixed 
points, Figs. [ 2 ] (b), (c) and (d). 

In Fig. [3] we show the mean field bifurcation diagrams 
of c given by Eq. (| 6 | as a function of the fraction of con¬ 
formists £ for e = 0.2 and e = 0. In this second case, 
there are absorbing states for large £ but for smaller val¬ 
ues the diagram hardly depends on the value of e. The 
leftmost vertical line at £ c marks the threshold at which 
the chaotic region ends. In Fig. [4] (a) we show an am¬ 
plication of the bifurcation diagram for small £ and in 
Fig. a (b) the corresponding Lyapunov exponent A and 
Boltzmann’s entropy 77 . For £ > £ c , A < 0 and 77 < 1/2. 
In Figs. [3] (a) and (b), the next vertical line at £ Q corre¬ 
sponds to the period-doubling bifurcation from a period 
two orbit to a fixed point. For £ c < £ < £ 0 the rather 
large fraction of contrarians causes symmetric oscillations 
of c. The rightmost vertical line at £{, corresponds to 
a pitchfork bifurcation from c = 1/2 to c > 1/2 when 


In Fig. [ 5 ] we show that Eq. ([ 7 ]) agrees with the numerical 
results for large connectivities k. 

In Fig. [6j we show bifurcation diagrams of the average 
opinion c as | J\ changes for small values of £. For £ = 0, 
Fig. [6] (a), there are no conformists and we have the case 
studied in Ref. [Tj. The bifurcation diagram is stretched 
horizontally as £ grows. For small \J\, c = 1/2, and as 
this parameter grows, there is a first bifurcation to a pe¬ 
riod 2 orbit. What seems like a second bifurcation to a 
period-four orbit, corresponds to two period-two bifurca¬ 
tions that depend on the initial opinion cq. Starting from 
below, the first and third branches (in red) correspond 
to co = 0.9, and the other two (in blue) to cq = 0.1. 


IV. SMALL-WORLD NETWORKS 

Real societies are not random, nor regular. It is inter¬ 
esting to study what happens when the topology changes, 
due for instance to advances in the transportation sys¬ 
tem, or to politics favoring mixing, etc. We studied the 
effect of rewiring a fraction p of links in a regular, one¬ 
dimensional society with connectivity k. This leads to 
the small-world networks first discussed by Watts and 
Strogatz m- As p , the long range connection probabil¬ 
ity, grows, small-world networks approach a mean field 
behavior. As we show in Fig. [7] the bifurcation diagrams 
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FIG. 4: (Color online.) (a) The mean field average opinion 
c as a function of the fraction of conformists £. For £ slightly 
smaller and also larger than £ c , we see four branches. Starting 
from below, the first and third branches correspond to Co = 
0.9 (in red), and the second and fourth to Co = 0.1 (in blue). 
After a 1,000 time steps transient, 32 consecutive iterations 
are plotted for each value of £ starting with Co = 0.9 and 
Co = 0.1. (b) The Lyapunov exponent A, bottom curve for 
small £ (in blue), and the entropy rj, top curve for small £ (in 
red), of the mean field average opinion c, Eq. (|6|, as functions 
of £. For each value of £, A is evaluated during 1,000 time steps 
and for r /, the unit interval is divided in 2 10 = 1, 024 equal size 
subintervals and the frequency with which each subinterval is 
visited is found during 100 x 2 10 = 102, 400 time steps after 
a 300 time steps transient. 



FIG. 5: (Color online.) Mean field critical values £ a as 

functions of k for | J\ = 5, circles (in red), and | J\ = 10 squares 
(in blue). The continuous curves are the graphs of Eq. the 
bottom one (in red) is for | J| = 5, and the top one (in blue) 
corresponds to | J\ = 10. The horizontal dashed lines show the 
asymptotic values (derived in the Appendix), £ a (5,oo) = 0.4, 
bottom dashed line (in red), and £ a (10, oo) = 0.45 top dahed 
line (in blue). Numerical data obtained after a transient of 
1,000 time steps starting with Co = 0.9. 


of c as functions of p are similar to those obtained by 
varying J in the mean field approximation. 

In Fig. [?] (a) we study a society of reasonable con¬ 
trarians. For small values of p , c fluctuates around 
c= 0.5. For slightly larger values, there are noisy oscilla¬ 
tions around two symmetric values, in a way reminiscent 
of the period-doubling bifurcations of deterministic sys¬ 
tems. For even larger values, p ~ 0.4, we observe the 
appearence of two different noisy oscillating states of pe- 



\J\ 


(b) 
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FIG. 6: (Color online.) Mean field bifurcation diagrams of 
the average opinion c as functions of | J| for several values of 
the fraction of conformists £. (a) £ = 0.0. (b) £ = 0.04. (c) 
£ = 0.07. (d) £ = 0.08. Starting from small |J| there is a 
first bifurcation to a period two orbit. For larger values of | J| 
there are four branches that lead. Starting from below, the 
first and third branches are obtained with an initial average 
opinion Co = 0.9 (in red), and the second and fourth with 
Co = 0.1 (in blue). For each value of |J|, after a transient 
of 1,000 time steps the next 32 iterations are plotted for two 
initial average opinions Co- 


riod 2. For Cq = 0.1 the orbit oscillates between the 
first and third branches, starting from the bottom (blue 
points), and for cq = 0.9 the orbit oscillates between the 
second and fourth branches (red points). This roughly 
corresponds to what is shown in Fig. [6] (a), although in 
that case the diagram is drawn as a function of | J| and 
here as a function of p. For £ = 0.05, Fig. [?] (b) we have 
a similar scenario and for larger values of £, Figs. 0(c) 
and (d), we find a great similarity with the mean field 
behavior. 

For £ = 0.2, Fig. [7] (c), c fluctuates around c = 0.5 
for small values of p and fluctuates around two branches 
for larger values. These branches agree with the period 2 
orbit of the mean field approximation, Eq. ([6]) shown in 
Fig. [3j For £ = 0.8, Fig. [7] (d), and p > 0.2, c fluctuates 
around one of two values, depending on cq , the average 
opinion at t = 0. These values also agree with the mean 
field approximation shown in Fig. [3] 

It is possible to roughly understand these results as¬ 
suming that the main contributions to the mean field 
character of the collective behavior come from the frac¬ 
tion of links that are rewired (long-range connections) 
that depend on p. The actual value of the field (2 h — 1) 
in Eq. Q is multiplied by a factor p , so that changing p 
is roughly equivalent to changing J. 

In Fig. [8] we show bifurcation diagrams of c on small- 
world networks as a function of the fraction of con¬ 
formists £ for some values of the long range connection 
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FIG. 7: (Color online.) Bifurcation diagrams of the average 
opinion c on a small-world network as a function of the long 
range connection probability p for different values of the frac¬ 
tion of conformists / with IV = 20, 000, | J\ =5, k = 20. The 
initial opinion of each agent is chosen at random between 0 
and 1 with Co = 0.1 (blue points) and Co = 0.9 (red points), 
(a) £ = 0.0. (b) £ = 0.05. (c) £ = 0.20. (d) £ = 0.8. In (c) 
and (d) the horizontal lines are the values of c of the period 
two orbit of Eq. <§• For each value of p , 32 points are plotted 
after a 1,000 time step transient for each Co- 


FIG. 8 : (Color online.) Bifurcation diagrams of the aver¬ 
age opinion c on a small-world network as functions of the 
fraction of conformists £ for different values of the long range 
connection probability p with N = 50,000, |J| = 5, k — 20. 
(a) p = 0.0. (b) p = 0.05. (c) p = 0.20. (d) p = 0.8. For 
each value of p, at t = 0 the opinion of each agent is chosen 
at random in such a way that the average opinion is Co = 0.1 
or co = 0.9. After a 50,0000 time step transient, 32 points 
are plotted. For 0.6 < £ < 0.8, the orbits that started with 
Co = 0.99 have c > 1/2, (points in red), and those that started 
with co = 0.1 have c < 1/2, (points in blue). 


probability p. This sequence of plots illustrate an un¬ 
expected behavior. For p = 0, Fig. [8] (a), there is a 
pitchfork bifurcation at £ ~ 0.9. For larger values of 
p, the pitchfork bifurcation occurs at smaller values of 
£ and an oscillating “bubble” is formed, Fig. [8] (b) for 
0.12 < £ < 0.22. This oscillating region grows with c and 
for p = 0.8, Fig. [8] (d), the bifurcation diagram is similar 
to the mean field one, Fig. 1(a)- 

We can explain this behavior assuming that the con¬ 
formist behavior promotes synchronization, as does p. So 
the system progressively synchronizes starting from high 
values of £, but this synchronization is not visible if the 
dynamics leads to fixed points. When the synchroniza¬ 
tion reaches the oscillating phases, it becomes manifest 
by means of the coherent oscillation of the population. 
So, this coherent dynamical behavior appears to start 
first in the vicinity of the bifurcation for £ ~ 0.2 and 
then, by increasing p, it propagates to lower values of £. 


V. SCALE-FREE NETWORKS 

In this section we present results of the model on un¬ 
correlated scale-free networks m- Starting from a fully 
connected group of m agents, other N — m agents join 
sequentially, each one choosing m neighbors among those 
already in the group. The choice is preferential, the prob¬ 
ability that a new member chooses agent i is proportional 
to its connectivity hi, the number of neighboring agents 


that agent i already has. Another way of building the 
network is choosing a random edge of a random node 
and connecting to the other end of the edge, since such 
an edge arrives to a vertex with probability proportional 
to kp(k) with p{k) the probability that a randomly se¬ 
lected node has connectivity k na. 

In Ref. [1] we showed that the dynamics of a model of 
a society whose agents are all reasonable contrarians on 
a scale-free network with m initially connected agents is 
comparable to the mean field approximation of Sec. Ill 
with connectivity k provided that 


k = i 


( 8 ) 


with a ~ 1.7 for scale-free networks with p(k) oc k 3 . In 
Fig- 0 we show that this result also holds for the model 
of societies studied here. In this Figure we compare the 
bifurcation diagrams, as £ changes, of the dynamics on 
a scale-free network with the corresponding mean field 
one, Figs. 9] (a) and (b), and the scale-free entropies with 
those of the mean field map, Figs. 9] (c) and (d). There 
is a reasonable agreement in both bifurcation diagrams. 
Both entropies show a good agreement where there is 
disorder, that is p > 1/2 but not when r) < 1/2. This can 
be understood from the bifurcation diagrams. While the 
scale-free network dynamics is stochastic and therefore 
the orbit visits many subintervals, the mean field one 
visits a smaller number. For example, for £ = 0.2, the 
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mean field dynamics has period two so 77 = 2/8 both in 
Figs. [9] (c) and (d). 



FIG. 9: (a) and (b). Bifurcation diagrams of the average 

opinion c as functions of the fraction of conformists £ on a 
scale-free network, large points (in blue) with m = 20 in (a) 
and m = 30 in (b), and the mean field approximation, Eq. (Im 
with k = 1.7 x 20 = 34 in (a) and k = 1.7 x 30 = 51 in (b), 
smaller points (in red), (c) and (d) Boltzmann’s entropy y 
of c as functions of ( on a scale-free network, top curve for 
f = 1/2 (in red), and of Eq. aGp, bottom curve for £ = 1/2 (in 
blue), (c) m = 20, k = 34. (d) m = 30, k = 51. In (a) and 
(b), for c on a scale-free network and for each value of £, at 
t = 0 the opinion of each agent is chosen at random in such a 
way that the average opinion is Co = 0.9 and for the mean field 
bifurcation diagram Co = 9 and Co = 0.1. After a transient of 
300 time steps, the next 32 values of c are plotted. For the 
entropy 77 and each value of £, the unit interval is divided in 
2 s = 256 equal size subintervals and the frequency with which 
each subinterval is visited is found during 100 x 2 8 = 25, 600 
time steps after a 300 time steps transient. 


VI. CONCLUSIONS 


In small-world networks with small values of p, Fig. [ 8 ] 
(b), the coherent oscillations appear first for a population 
with a small fraction of conformist rather than for a pure- 
contrarian one. This is probably due to the fact that 
coherent oscillations are a signal of synchronization, and 
the presence of conformists increases the synchronization. 
On the other hand, if the synchronized dynamics leads 
to a stable fixed point, the degree of synchronization is 
not manifest. As a result, the first synchronized zone is 
near the bifurcation point £ ~ 0 .. 

Dynamics on small-world networks approaches that of 
the mean field approximation as the long range connec¬ 
tion probability p grows. Dynamics of c on scale-free net¬ 
works is similar to that of the mean fileld approximation 
provided k = am with a ~ 1.7. 
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Appendix 

The bifurcation points £ 0 and £ b can be found from the 
mean field evolution, Eq. © , since at both values, the 
absolute value of the derivative of this expression must 
be 1. For large k the mean field approximation for the 
average opinion, can be approximated by |T 


The dynamics of the mean field approximation of the 
average opinion c, Eq, © as a function of £ is chaotic 
with periodic windows when 0 < £ < £ c and oscil¬ 
lates periodically between two symmetric values when 
Cc < ? < C,- For £ a < £ < &>, c = 1/2, and for < £, 
c > 1/2 (c < 1/2) if c 0 = 0.9 (c 0 = 0.1). Also £ a ~ 1 - & 
and £ a approaches a limit value value as k grows. Eq, © 
depends on the connectivity k and on the transition prob¬ 
ability r that in turn depends on k, \ J\, e, and q. As far as 
we have explored the parameter space, the above descrip¬ 
tion is generic with the exception that for large values of 
k the chaotic phase is bounded below by ^,1 > 0. For 
0 < £ < there are two period two orbits that depend 
on Co and for ^ < £ < £ c the orbits are chaotic with 
periodic windows. 


d = I dx\ 


27tc(1 — c) 


exp 


k(x — c ) 2 
2 c(l — c) 


[£t(z; J) + {1-£)t(x;-J)} . 


(9) 


Expanding the right-hand side term of this expression 
around c = 1/2 up to first order 

d = y— J dxex p (—2 k(x — c) 2 ) • 

[£t(x; J) + (l-£) t(x-,-J)} 

= \J^~ Jdyexp(-2ky 2 )- 

[C T (y + c; J) + (1 — ^)r(y + c; — J)]. (10) 
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By approximating 


We denote Eq. (10 1 by 

To proceed we need the derivative of r given by Eq. Q 
near c = 1/2 and small y. From Eq. Q 


f(y; J) ~ J(1 - 4JV) ~ J exp(-4 J 2 y 2 ) 




4 J exp(—4 Jy) 


v—i/2+y (l + exp(-4 Jy)) 2 ' 
We note that f(y ; — J) = —f(y; J )• Then 

g{^J)= / dy\j^exv(-2ky 2 ) ( 2 £ - l)/(y; J). 


we get for £ a 


(2& - 1) J 


jfe 


2 J 2 + k 


= -1 


i.e., 


It is now straightforward to check that 


5 (1 - 6 J)= fl (£; -J) = -<?(£; J). 





If 3 (^ a ; J) = —1 then 3 (^ 5 ) = 1 with ^ = 1 — £ a and the 

the two bifurcation points are symmetric with respect to The limit of this last expression when k —> 00 is 

£ = 1/2. £»(>/ 00 )) = ( 1 / 2)(1 — 1 /J). 
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